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Abstract 

Branching Processes in Random Environment (BPREs) (Z n : n > 0) are the generalization 
of Galton- Watson processes where in each generation the reproduction law is picked randomly 
in an i.i.d. manner. In the supercritical case, the process survives with positive probability and 
then almost surely grows geometrically. This paper focuses on rare events when the process 
takes positive but small values for large times. 

We describe the asymptotic behavior of P(l < Z n < k\Zo = i), k, i G N as n — > oo. More 
precisely, we characterize the exponential decrease of ¥(Z n = k\Zo = i) using a spine repre- 
sentation due to Geiger. We then provide some bounds for this rate of decrease. 
If the reproduction laws are linear fractional, this rate becomes more explicit and two regimes 
appear. Moreover, we show that these regimes affect the asymptotic behavior of the most 
recent common ancestor, when the population is conditioned to be small but positive for large 
times. 

AMS 2000 Subject Classification. 60J80, 60K37, 60J05, 60F17, 92D25 

Key words and phrases, supercritical branching processes, random environment, large deviations, 
phase transitions 

1 Introduction 

A branching process in random environment (BPRE) is a discrete time and discrete size population 
model going back to |33l 17]. In each generation, an offspring distribution is picked at random, 
independently from one generation to the other. We can think of a population of plants having 
a one-year life-cycle. In each year, the outer conditions vary in a random fashion. Given these 
conditions, all individuals reproduce independently according to the same mechanism. Thus, it 
satisfies both the Markov and branching properties. 

Recently, the problems of rare events and large deviations have attracted attention [2"5l 191 H2l |2"91 |2"51 
[TU] , However, the problem of small positive values has not been treated except in the easiest case 
which assumes non-extinction, i.e. P(Zi > 1\Zq = 1) > 1 (see [9]). For Galton- Watson processes, 
the explicit equivalent of this probability is well-known (see e.g. [8] [Chapter I, Section 11, Theorem 



3]). In particular, denoting by / the probability generating function of the reproduction law, we 
have for k large enough, 



g := lim i logP(l < Z n < k\Z = 1) = /'(p e ), where Pe = ¥(3n e N : Z n = | Z = 1). (1.1) 

n— >-oo 

Moreover, the rate of decrease remains equal to g if k n decreases sub-exponentially. It means that 
as soon as k n / exp(Sn) — > as n — > oo for every S > 0, then lim n _>. 00 i logP(l < Z n < k n \ Zq = 
1) = g. In this paper, we focus on the existence and characterization of g in the random environ- 
ment framework. It is organized as follows. 

First, we give the classical notations and properties of BPRE. In the next section, we state our 
results. We prove the existence of g and a characterization of its value via a spine construction, 
give a lower bound and an upper bound which have natural interpretations. Finally, we specify 
our results in the linear fractional case, where two regimes appear, which are also visible in the 
time of the most recent common ancestor (MRCA). 

In the rest of the paper, the proofs of these results are presented. Section [3] deals with a tree 
construction due to Geiger, which is used in Section [4] to characterize p. In Section 5.2 we prove 



that p > under suitable assumptions. In Section |5.3[ we prove a lower bound for g in terms of 
the rate function of the associated random walk. Finally, in Sections |6.1| and [6. 2| the statements for 
the linear fractional case are proved using the general results obtained before, whereas in Section 
[7J we present some details on two examples. 

For the formal definition of a branching process Z in random environment, let Q be a random 
variable taking values in A, the space of all probability measures on N = {0, 1,2,.. .}. We denote 

by 

m q :=J2k q({k}) 

k>0 

the mean number of offsprings of q 6 A. For simplicity of notation, we will shorten q({-}) to 
q(-) throughout this paper. An infinite sequence £ = (Qi, Q%, ■ ■ .) of independent, identically 
distributed (i.i.d.) copies of Q is called a random environment. Then the integer valued process 
(Z n : n > 0) is called a branching process in the random environment £ if Zq is independent of £ 
and it satisfies 

C(Z n | £, Z ,..., Z„_i) = Q* n z - 1 a.s. (1.2) 

for every n > 1, where q* z is the z-fold convolution of the measure q. We introduce the probability 
generating function (p.g.f) of Q n , which is denoted by /„ and defined by 

oo 

Ms) :=5> fc Q„(fc), (sG[0,l]). 

fc=0 

In the whole paper, we denote indifferently the associated random environment by £ — (/i, /2, ...) 
and £ = (Qi, Q2, •••)• The characterization (|1.2| of the law of Z becomes 



E[s z "\£, Z ,...,Z n _ 1 ]=Ms) Zn - 1 a.s. (0<«<l,n>l). 

Many properties of Z are mainly determined by the random walk associated with the 
environment (S n : n g No) which is defined by 

5*0 = 0, S n — S n -i=X n (n > 1), 
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where 

X n := \ogm Qn = log/4(l) 

are i.i.d. copies of the logarithm of the mean number of offsprings X := log(mg) = log(/'(l)). 
Thus, one can check easily that 

E[Z n \Q 1 ,...,Q n ,Z = l] = e s " a.s. (1.3) 

There is a well-known classification of BPRE [7], which we recall here in the case E[|X|] < oo. In 
the subcritical case (E[X] < 0), the population becomes extinct at an exponential rate in almost 
every environment. Also in the critical case (E[X] = 0), the process becomes extinct a.s. if we 
exclude the degenerated case when F(Z\ = 1\Zq = 1) = 1. In the supercritical case (E[X] > 0), 
the process survives with positive probability under quite general assumptions on the offspring 
distributions (see 33 ). Then E[Zi log + (2'i)//'(l)] < oo ensures that the martingale e~ Sn Z n has 
a positive finite limit on the non-extinction event: 

lim er s "Z n = W a.s., F(W > 0) = P(Vn e N : Z n > 0\Z Q = 1) > 0. 

n— s-oo 

The problem of small positive values is linked to the left tail of W and the existence of harmonic 
moments. In the Galton- Watson case, we refer to [51 E21 [XZ] ■ For BPRE, Hambly [22] gives the 
tail of W in 0, whereas Huang & Liu [23, 24 have obtained other various results in this direction. 



2 Probability of staying bounded without extinction 

Given the initial population size k, the associated probability is denoted by Pfc(-) := P(-|Zo = 
For convenience, we write P(-) when the size of the initial population is taken equal to 1 or does 
not matter. Let fi t „ be the probability generating function of Z n started in generation i < n : 

fi,n ■= fi+1 ° fi+2 ° ■ ■ ■ ° fn> fn.n = Id a.S. 

We will now specify the asymptotic behavior of ¥i(Z n = j) for i,j > 1, which may depend both on 
i and j. One can first observe that some integers j cannot be reached by Z starting from i owing 
to the support of the offspring distribution. 

The first result below introduces the rate of decrease q of Fi(Z n = j) for i, j > 1 and gives a 
trajectorial interpretation of the associated rare event {Z n = j}. The second one provides general 
conditions to ensure that q > 0. It also gives an upper bound of g, which may be reached, in 
terms of the rate function of the random walk S. This bound corresponds to the environmental 
stochasticity, which means that the rare event {Z n = j} is explained by rare environments. The 
next result yields an explicit expression of the rate g in the case of linear fractional offspring dis- 
tributions, where two supercritical regimes appear. The last corollary considers the most recent 
common ancestor, where a third regime appears which is located at the borderline of the phase 
transition. 

Let us define 

!:={]> I : P(Q(j) >0,Q(0) > 0) > 0} 
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and introduce the set Cl({z}) of integers that can be reached from z £ I, i.e. 

Cl({z}) := {k > 1 : 3n > with F z (Z n = k) > 0}. 

Analogously, we define Cl(I) as the set of integers which can be reached from I by the process Z. 
More precisely, 

Cl(l) := {k > 1 : 3n > and j £ 1 with Pj(Z n = k) > 0}. 

We observe that 1 C Cl{T) and if P(Q(0) + Q(l) < 1) > and P(Q(0) > 0, Q(l) > 0) > 0, then 
CZ(I) = N. 

We are interested in the event {Z„ = j} for large n. First, we recall that the case f{Z\ = 0) = 
is easier and the rate of decrease of the probability is known [5] . Indeed, then Z is nondecreasing 
and for k > j £ N such that Vk(Zi — j) > for some I > 0, we have 

lim ±\ogP k (Z n =j) = -kg, with e = -logPi(Zi = 1). 

n— »oo 

We note that in the case P\(Zi = 1) = P(Zi = 0) = 0, the branching process grows exponentially 
in almost every environment and the probability on the left-hand side is zero if n is large enough. 
Thus, let us now focus on the supercritical case with possible extinction, which ensures that T is 
not empty. The expression of p in the next theorem will be used to get the other forthcoming 
results. 

Theorem 2.1. We assume that E[X] > and V{Z\ = 0) > 0. Then the following limits exist and 
coincide for all k,jE Cl{I), 

g := - lim ±lo E P k (Z n =j) = - lim I logE[Q n (z )/o,«(0) Zn - 1 n^ 1 1 /( (&„(<)))] 

where z$ is the smallest element of I. The common limit q belongs to [0, oo). 

The proof is given in Section [4] and results from Lemmas |4.1| and |4.2| 

The right-hand side expression of g correspond to the event {Z n = j} explained by a "spine struc- 
ture" . More precisely, one individual survives until generation n and gives birth to the j survivors 
in the very last generations, whereas the other subtrees become extinct (see forthcoming Lemma 



3.2 for details). However, we are seeing in the linear fractional case (Corollary 2.3) that a multi- 



spine structure can also explain {Z n — j} in some regime. Thus the optimal strategy is nontrivial 
and will here only be discussed in the linear fractional case. 



The proof of Theorem 2.1 is easy if we consider the limit of ^ logPi(Z„ = 1) as n —> oo. In this 
case, a direct calculation of the first derivative of /o jTl yields the claim. However, the proof for 
the general case is more involved. Here, we use probabilistic arguments, which rely on a spine 
decomposition of the conditioned branching tree via Geiger construction. 



We also note that we need to focus on i,j £ Cl(T). Indeed, limn.^ ^ logPi(Z„ = i) and 
lim n _>. 00 - logPi(Z n = j) may both exist and be finite for i ^ j, but have different values. To see 
that, one can consider two environments q\ and qi such that 

P(Qi = <Zi) = 1 - P(Qi = 92) > 0; gi(l) = l; ga(0) + <a(2) = l. 

Moreover, the case —00 < limsup„_ j . 00 - logP fc (Z„ = k) < lhninf„_j. 00 ~ log Pi (Z n = k) < with 
k > 1 is also possible. These results are developed in the two examples given in Section [7] at the 
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end of this paper. 

In the Galton- Watson case, / is constant and for every i > 0, = f a.s. Then /i, n (0) — > p e as 



oo and we recover the classical result (1.1) 



The results and remarks above could lead to the conjecture g = — logE[/'(p(/))], where p(f) = 
inf{s £ [0, 1] : f(s) = s}. Roughly speaking, it would correspond to integrate the value obtained 
in the Galton- Watson case with respect to the environment. The two following results show that 
this is not true in general. 

To prove that the probability of staying small but alive decays exponentially (i.e. p > 0) requires 
some assumptions. To avoid too much technicalities, we are assuming 

Assumption 1. There exists 7 > such that Q(0) < 1 — 7 a.s. and E[|X|] < 00. 

Similarly to give an upper bound of g in terms of the rate function of the random walk S, we 
require the following Assumption. The non-lattice condition is only required for a functional limit 
result which is taken from [5], whereas the truncated moment assumption is classically used for 
lower bounds of the survival probability of BPREs. 

Assumption 2. We assume that S is non-lattice, i.e. for every r > 0, P(X £ rZ) < 1. 

Moreover, there exist e > and a £ N such that for every x > 0, 

E[(log+ t Q (a)) 2+£ \X > -x] <0 o, 
where log + x := log(max(x, 1)) and £ q (a) is the truncated standardized second moment 

oo 

f,(o) :=J2y 2 l(y)/ m l a£N,q£A. 

y=a 

Proposition 2.2. We assume that there exists s > such that E[e~ sX ] < oo. 

(i) If Assumption^ is fulfilled, then p > 0. 

(ii) IfV(X > 0) = 1 or Assumption^ holds, then 

g < - log inf E[exp(AX)]. 
We note that the exponential moment assumption is equivalent to the existence of a proper rate 



function A for the lower deviations of S. The lower bound (i) is proved in Section 5.2 The second 
bound is the rate function of the random walk S evaluated in 0, say A(0). Indeed, recalling that 
E[X] > 0, the supremum in the Legendre transform can be taken over M + instead of M. Extract- 



ing — 1 yields the upper bound above. It is proved in Section 5.3 and used for the proof of the 
next Corollary |2.3[ It can be reached and has a natural interpretation in terms of environmental 
stochasticity. Indeed, one way to keep the population bounded but alive comes from a 'critical en- 
vironment', which means S„ ~ 0. Then E[Z„ | £] — exp(S* n ) is neither small nor large and one can 
expect that the population is positive but bounded. The event {S n ~ 0} is a large deviation event 
whose probability decreases exponentially with rate A(0). This bound is thus directly explained 
by the environmental stochasticity. 
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Now, we focus on the linear fractional case. We derive an explicit expression of g and describe 
the position of the most recent common ancestor of the population conditioned to be positive but 
small. We recall that a probability generating function of a random variable R is linear fractional 
(LF) if there exist positive real numbers m and b such that 

/OO = 1 - 



„i- 1 + 6m- 2 (l-s)/2' 

where m = /'(l) and b = /"(l). This family includes the probability generating function of 
geometric distributions, with b — 2m 2 . More precisely, LF distributions are geometric laws with a 
second free parameter b which allows to change the probability of the event {R = 0}. 

Corollary 2.3. We assume that f is a.s. linear fractional, E[\X\] < oo, E[X 2 e~ x ] < oo and 
¥{Zi = 0) > 0. 

We assume also that either P(A" > 0) = 1 or Assumption^hold. Then 



-logE[e- x ] , if E[Xe" x ]>0 

-loginf A > E[exp(AX)] (= A(0)) , else 



(2.1) 



This result is also stated in the PhD of one of the authors and can be found in [IT] . On the level 

of large deviation (log scale), two regimes in the supercritical case are visible. 

If E[Xe- x ] < 0, the event {Z n = k} is a typical event in a suitable exceptional environment, say 

'critical'. This rare event is then explained (only) by the environmental stochasticity. 

If E[ATe~ x ] > 0, we recover a term analogous to the Galton- Watson case, which is smaller than 

A(0). The rare event is then due to demographical stochasticity. 

These two regimes seem to be analogous to the two regimes in the subcritical case, which deal with 
the asymptotic behavior of Z n > 0, see e.g. [21 [21] [20] ■ Such regimes for supercritical branching 
processes have already been observed in |26j in the continuous framework (which essentially rep- 
resents linear fractional offspring-distributions). 

Let us now focus on the most recent common ancestor (MRCA) of the population conditionally on 
this rare event. More precisely, let T n be the set of all ordered, rooted trees of height exactly n. 
We refer to |31j for classical definitions. We say that an individual i„ 2 in generation n 2 > rii stems 
from an individual i ni iff there are individuals i n2 _i, . . . , i ni +i such that i ri2 is a child of i„ 2 _i, 
i n2 —i is a child of i n2 — 2i ■ ■ ■ an d i ni +i is a child of i ni . Let T n G T n be the random branching 
tree, generated by the process (Zk)o<k<n conditioned on Z n > and denote by MRCA n the most 
recent common ancestor of the population at time n. More precisely, we consider the events 

A]?: := { all individuals in generation n stem from one individual in generation n — k} 

and define the age of the MRCA in generation n as the number of generations one has to go back 
in the past until all individuals in generation n have a single common ancestor : 

MRCA n := min{fc = 1, 2 . . . , n \A% holds }. 

The case P(Zi = 0) = is trivial : extinction is not possible, so q = — logPi(Zi = 1) and 
P2 (MRCA n = n\Z n = 2) = 1. It is excluded in the statement below. Moreover, for the sake of 
simplicity, we restrict the study of the MRCA to starting from one individual and conditioning on 
Zn = 2. 
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Corollary 2.4. We make the same assumptions as in the previous corollary. 

(i) IfE[Xe- x ] < 0, then 

]xmm{Fi(MRCA n =n\Z n = 2) >0 ; liminf Pi (MRCA n = l\Z n = 2) > 0. 

(ii) IfK[Xe~ x ] = ? then for every sequence (x n )neN such that x n € /or every n 7 we have 

lim i logPi (Mi?CA„ = x„|Z„ = 2) = 0. 

n— »oo 

Moreover, under the additional assumption E[/"(l)] < oo, there exist two positive finite constants 
c, C such that for every n £ N, 

c < n Fx(MRCA n = n\Z n = 2) < C 

and if E[/"(l)/(l — /(0)) 2 ] < oo, £/ien /or every 8 £ (0, 1), £/iere exzsi iwo positive finite constants 
c', C" smc/i f/iai for every n £ N, 

c' < n 3/2 Pi(MECA» = [fc] |Z„ = 2) < C". 

(Hi) IfE{Xe~ x ] > anrf E[e^ 1_s ^] < oo for some s > 0, then for every S e (0, 1], 

lim sup i \og¥ x {MRCA n > 5n\Z n = 2) < 0. 

Thus three regimes appear for the most recent common ancestor of the population. 
If E [Xe _x ] > 0, which we call the 'strongly' supercritical case, the MRCA is at the end (close to 
the actual time). The probability that the MRCA is far away from the final generations decreases 
exponentially. Such a result is classical for branching processes which do not explode, such as 
subcritical Galton- Watson processes conditioned on survival. It corresponds to a spine decom- 
position of the population whose subtrees become extinct Conditionally on {Z n = 2}, S is 
still a random walk with positive drift and will be typically large. Thus the conditioned process 
is typically small throughout all generations (as in the Galton- Watson case) as growing and then 
becoming small again within the favorable environment has a very small probability. Consequently, 
the MRCA will be close to generation n. 

But in the 'weakly' supercritical case (E[Xe _x ] < 0), conditionally on {Z n — 2}, the MRCA is 
either at the beginning (close to the root of the tree) or at the end (close to generation n). Such 
a situation is much less usual. It has already been observed in [TH] for the subcritical reduced tree 
of linear fractional BPRE conditioned to survive. Here, as indicated in the proof of Proposition 
|2.2| (ii), the random walk S conditioned on {Z n = 2} typically looks like an excursion. It means 
that S is conditioned on the event {min{5o, . . . S n } > 0, S n < c}. In such an environment, sub- 
trees that are either born at the beginning or at the end may survive until the end. All subtrees 
being born at some generation |^<5rzJ , S £ (0, 1) experience an unfavorable environment and become 
extinct. This can be seen as follows. During an excursion from to n, typically S\g n \ 3> and 
thus e s "~ s L<s™j <§; and the corresponding subtree will become extinct. 

Finally, in the intermediate case (E[Ae _x ] = 0), the MRCA is close to the end, but the proba- 
bility that the MRCA is far away from the end only decreases polynomially. The intermediately 
supercritical regime is in-between the two regimes described above and conditioned on {Z n = 2}, 
the typical environment will neither be an excursion nor a random walk with positive drift. 
One can expect several more detailed results describing the three regimes, which are beyond the 
scope of this paper. 
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3 The Geiger construction for a branching process in vary- 
ing environment (BPVE) 



In this section, we work in a quenched environment, which means that we fix the environment 
e := (gi, (72, • • •)■ We consider a branching process in varying environment e and denote by P(-) 
(resp. E) the associated probability (resp. expectation), i.e. 

P(Zi = fci, • • • , Z n = k n ) = V(Zi = fci, • • • , Z n = k n \£ = e). 

Thus (/i, /2, . . .) is now fixed and the probability generating function of Z is given by 

E[s z " \ Z = k]= fo^ n (s) k (0<s<l). 

We use a construction of Z conditioned on survival, which is due to [18] [Proposition 2.1] and extends 
the spine construction of Galton- Watson processes |25j . In each generation, the individuals are 
labeled by the integers i = 1, 2, ••• in a breadth-first manner ('from the left to the right'). We 
follow then the 'ancestral line' of the leftmost individual having a descendant in generation n. This 
line is denoted by L. It means that in generation fc, the descendance of the individual labeled 
survives until time n, whereas all the individuals whose label is less than become extinct before 
time n. The Geiger construction ensures that to the left of L, independent subtrees conditioned 
on extinction in generation n are growing. To the right of L, independent (unconditioned) trees 
are evolving. Moreover the joint distribution of h k and the number of offsprings in generation k is 
known (see e.g. [1 for details, where L := L — 1) and for every k > 1, z > 1 and 1 < 11 < z, 

P(7 n 717 17N nl ( , P(Z n >0\Z k ^l)P(Z n ^O\Z k ^iy- 1 
P(Z k = z,U = l\Z k . x = 1, Z n > 0) = q k {z) P{Zn>0 \ Zk _ 1 = 1) ■ 



Let us give more details of this construction. We assume that the process starts with Zq — z and 
denote P z (-) := P(-\Z = z). We define for < k < n, 

Pk,n := P(Zn > | Z k = 1) = 1 - /fc,„(0), Pn,n ■= 1- 

We can specify the distribution of the number Y k of unconditioned trees founded by the ancestral 
line in generation k, at the right of L^. In generation 0, for < i < z — 1, 

P Z {Y Q = i\Z n > := P z L = z-i\ Z n >0) = — — r 

P(Z n > \Z = z) 

_ l-/o,n(0) 



P z {Z n > 0) 



hnW ■ (3-2) 



More generally, for all 1 < k < n and i > 0, (3.1) yields 

P{Y k = i\Z n > 0) := P{Z k - L k = i\Z n > 0, Z k ^ = 1) 

OO 

= P(z k =jX k =j-i\z n >o,z k _ 1 = i) 



j=i+l 

Pfe.n 



E zkWkAoy-*- 1 - (3-3) 



' .7=1+1 

Finally, we note that f n , n (0) — 0, thus for k = n, we have P(Y n = i\Z n > 0) = 2sl&±i)_ 
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Figure 1: Geiger construction with trees conditioned on extinction and unconditioned 
trees. 

Here, we do not require the full description of the conditioned tree since we are only interested in 
the population alive at time n. Thus we do not have to consider the trees conditioned on extinction, 
which grow to the left of L. We can construct the population alive in generation n using the i.i.d 



random variables Y , Yi, Y" 2 j • • • j Y-n whose distribution is specified by (3.2 1 and (3.3 1 : 

F(Y k = i) = P(Y k = i\Z n > 0) 

Let (Zj k ^)j>o be independent branching processes in varying environment which are distributed 
as Z for j > k and satisfy 

Zf := for j < k, Zf ] := Y k . 
More precisely, for all < k < n and zq, ■ ■ ■ ,z n >0, 

P(Z ( k) = 0, • • • , Z [ k k \ = 0, Z { k k) = z k ,Z ( k % = z k+1 , z n ) 

= P{Y k = z k )P(Z k+1 = z k+ i, ■ ■ ■ , Z n = z n | Z k = z k ). 

The sizes of the independent subtrees generated by the ancestral line in generation k, which may 
survive until generation n, are given by (Z^)o<:j<: n , < k < n — 1. In particular, 

C(Z n \Z n > 0) = £(ZW + ... + Z^- 1] +Yn + 1). (3.4) 

Lemma 3.1. The probability that all subtrees emerging before generation n become extinct before 
generation n is given for z > 1 by 

n— 1 n— 1 

p,(i(°) + . . . + z<?-v = o) = [] P*(z ( n k) = o) = ^ n #(A.n(Q)). 

P-1 n 

fc=0 r ' k=0 
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where we use the following convenient notation fo(s) := s z , P-i, n ■— Pz(Z n > 0). 

Proof. First, we compute the probability that the subtree generated by the ancestral line in gen- 
eration k does not survive until generation n, i.e. {Zn — 0}. By (3.3), for k > 1, 

oo 

P(^ fe) =0) - ^P z (f fe -z)P(Z„-0|Z fc = l ) 

i=0 

oo oo 



Pk.n 



E E ft(j)/fc,n(0) J - i - 1 • 



OO oo 

^E E ft(j)/fc,n(0y- 1 

Pfe— 1,71 - r, - -ii 



OO 

Pfc,n 



E-^W-M )'" 1 



P& l 5 'f 7 = 1 

-^-/fc(/*,n(0)). (3-5) 

Pfc-l,n 



Similarly, we get from (3.2) that 



p z (z^ = o) = Y,p z (Y = i\z n >o)p(z n = o\z 



t=0 

= ^Wor' = ^=-/S(A,„(o)) 

P-l,n P-l,n 

with the convention /o( s ) = s? ' '• Adding that the subtrees given by {Zj k )j>o are independent, a 
telescope argument yields the claim. □ 

For the next lemma, we introduce the last generation before n when the environment allows 
extinction : 

K n := sup{l < k < n : q k (0) > 0}, (sup0 = O). 
Note that K n only depends on the environment up to generation n. 
Lemma 3.2. Let zq be the smallest element in X. Then, 

p zo (z n = zo i z n > o) = x n ^^-f k (h^M) x n 

P t\ / / - 1 .Kvi i i-i P^ 1 . K, ri ii 

K—V ' J — ^n + i- 

where we recall the following convenient notation fo(s) = s z , P-i, n — P z(Z n > 0). 

Proof. By definition of X and Zq, q(0) > implies q{k) = for every 1 < k < z n . We first deal 
with the case n n > 0. Then, 

9«„(0)>0, q Kn (k) = 0ii l<k<z ; <W+i(0) = • • • = g„(0) = 0. 

In particular the number of individuals in generation n n is at least zq times the number of individ- 
uals in generation n n — 1 who leave at least one offspring in generation n n . Moreover, as extinction 
is not possible after generation K n , it holds that Z Kn < Z Kn+ i < • • • < Z n . 

Let us consider the event Z n — z > 0. Then Z K \ > and Z Kn > z . So Z Kn = Z Kn+ i = 
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• • • = Z n = zq and only a single individual in generation n n — 1 leaves one offspring (or more) 
in generation k„. This individual lives on the ancestral line. Thus all the subtrees to the right 
of the ancestral line which are born before generation n n have become extinct before generation 
k„, i.e. ZkJ = . . . = zj^ n 1 ' = 0. In generation K n — 1, the individual on the ancestral line has 
zq offsprings and Y Kn = zo — 1. After generation «„, all the individuals must leave exactly one 
offspring to keep the population constant until generation n, since g K „+i(0) = • • • = 9n(0) = 0. 
This probability is then given by qj(l) z ° in generation j > n n . Moreover, (3.3) simplifies to 
P(Y Kn = zq — 1) = (7 Kii (zo)/Pk„-i,k„- Using the previous lemma, it can be written as follows: 




Figure 2: Illustration of the proof of Lemma [3~2 



P z {Z n — Zq) 

= p Zo (zjv = ... = z^ = o)P(r K 



l)P Z0 {Z^ + ■■■ + Zlr 1] +Yn + 1 = Zo) 



n 



Pfc.n 



fe=0 

g«„(^o) 

Pk„-1,k„ 



Pfc-l,n 

K.-l 



/*(/*,«„ («)) 



P «.„-!,«„ 



n 

fe=0 



Pfc-l,ri 



#(/m»(«)) 



P Z0 + . . . + Zi"- 1 ' + Y n + 1 = z ) 

n 



Recall that after generation n n , each individual has at least one offspring and thus pj <n = Pj >Kn for 
any j < K n . This ends up the proof in the case K n > 0. The case when n n = is easier. Indeed, 



since g K „+i(0) 



P Z0 (Z n = z ) = P Z0 (Z 1 = --- = Z n = Zo) = q 3 (l) Z0 

3=1 

- 9n(0) = and Z is nondecreasing until generation n. 



□ 



4 Proof of Theorem |2.1| : the probability of staying positive 
but bounded 



In this section, we prove Theorem |2.1| with the help of two lemmas. The first lemma establishes 
the existence of a proper 'common' limit. 

Lemma 4.1. Assume that z > 1 satisfies P(Q(0) > 0, Q(z) > 0) > 0. 
Then for all k,j £ Cl({z}), the following limits exist in [0, oo) and coincide 



lim ilogP fc (^ n = j) = lim ifogP z (Z„ = z). 
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Moreover, for every sequence k n such that k n > z for n large enough and k n /n — > as n — > oo, 

lim \ logP,(Z„ = z)= lim i logP z (l < Z„ < fc„). 

Proof. Note that for every fc > 1, Pfe(^i = z) > since 

P fc (Z! = « | Qi) > gi(0) fe - 1 Q 1 (z), P(Q(0) > 0, Q(z) > 0) > 0. 

We know that by Markov property, for all m, n > 1, 

F z (Z n+m = z)> ¥ z (Z n = z)P z (Z m = z). (4.1) 

Adding that F z (Zi — z) > 0, we obtain that the sequence (a n ) ne ^ defined by a n :— — logPz(Z„ = 
z) is finite and subadditive. Then Fekete's lemma ensures that lim n _ > . 00 a n /n exists and belongs to 
[0, oo). Next, if j, k e Cl({z}), there exist I, m > such that F z {Z t = j) > and P z (Z m = fc) > 0. 
We get 

F k {Z n+l+l = j) > P fc (Zi - z)P z (Z„ = z)P,(Z ; = j) 

and 

F z (Z m+n+1 =z)> F z (Z m = k)P k (Z n = j)P J (Z 1 = z). 
Adding that Fj{Z\ = z) > 0, we obtain 

liminf i logP fe (Z„ = j) > lim ± logP z (Z n = z) > lim sup I logP fe (Z„ = j), 
which yields the first result. 

For the second part of the lemma, we first observe that F z (Z n = z) < P z (l < Z n < k n ) for n large 
enough. To prove the converse inequality, we define for e > the set 

A £ := {q € A|g(0) > e, q(z) > e}. 

According to the definition of I and the assumption, F(Q G A £ ) > if e is chosen small enough. 
Thus, we get 

F z (Z n = z) > F z (l < Z n _! < fc„) mm F^Z, = z) 

l<j<k n 

>P Z {1< Z n ^<k n )P(Q eA e ) min ^ l {Z l =z)P 1 {Z 1 = Q\Qy- 1 \Q&A e ] 

l<j<k n 

>F Z (1< Z„_x < fc„)P(Q g A) £ fc ". 
Using the logarithm of this expression and letting n — > oo yields 

lim i logP z (Z„ = *) > limsup (i logP z (l < Z n _i < fc„) + log(e)^) . 
Adding that k n = o(n) by assumption gives the claim. □ 

Now, we prove a representation of the limit p in terms of generating functions. It will be useful in 
the rest of the paper. 

Lemma 4.2. Assume that Pi(Zi = 0) > 0. Then for all i,j e Cl(l), 

lim MogP j (Z„=j)= lim MogEfQ„(^o)/o,„(0) 2o - i n/i(/^(°)) 

i=l 

where zq is the smallest element in X. 
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We note that Pi(i>i = 0) > is equivalent to P(Q(0) > 0) > and in view of Lemma 4.1 we only 
have to prove the result for k — -- j — - zq, where zq is the smallest element in X. Differentiation of 
the probability generating function of Z n yields the result for Zq = 1. The generalization of the 
result for zq ^ 1 via higher order derivatives of generating functions appears to be complicated. 
Instead, we use probabilistic arguments, involving the Geiger construction of the previous section. 

Proof. First, the result is obvious when z = 1 e 1 since 

d n ~ ^ 

P{Z n = 1\£) = -/ 0) „( s )| s=o = /;( ) • n fKfUO)). 



i=l 



For the case zq > 1, we start by proving the lower bound. Using (3.4 1, Lemma 3.1 (3.2) and 
recalling that P Zg (Z n > 0\£) — p-i >n , we have 

F Z0 (Z n = z ) = e[p 2q (Z„ = z \Z n > 0,£)P z „{Z n > 0\£) 

= E\F Za (Z^ + ■■■ + Z^-V +Y n + 1 = z a \£)P Zo (Z n > 0\£) 
I P Z0 (Z^ + ■■■ + Z^ = 0, Y n = z - l\£)P Zo (Z n > 0\£) 



E 



= E 



,(Zn>0\£) 



Pn—l,n 
Pn-1,« P-l. 



E[/<(An(0)) 



i=0 



n-1 



QMUfiifiAv) 

i=0 

Recalling also that / (s) = zos 2 " -1 , we get 

n 

hm i logP Z0 (Z„ = z ) > limsup i fogE^o)/^ 1 ^) [] /f(/i,„(0)) 



(4.2) 



i=l 



Let us now prove the converse inequality. Following the previous section, zq is the smallest element 
in X and k„ is the (now random) last moment when an environment satifies Q(0) > 0. We 
decompose the event {Z n = zq} according to K n : 



n 

\(Z n = zq) = ^e[p 2[) (Z„ = z \£,Z n > 0)P Zo {Z n > 0\£);K n 

fc=0 



= k 



Using that conditionally on n n — k, P Zo (Z n > 0\£) = P Zo (Z k > 0\£) and Lemma 3.2 we get by 
independence 



v zo (z n = zo) =j2^o(z k > q\£) {zo) Pfc " i,fc n fKfiMo)) n ; k - 

to L P^ 1 ^ p -^ fc to j=k+l 

n k— 1 n 

<^E[Q fe (z )n/i(/a(o))] n ^[QiW ] 

k=0 i=0 j=k+l 

n k— 1 

= $>[Q fc (2 ) I] /i(/i,*(0))]E[Q(l)»°] n - fc - 1 +E[Q(iy 



k 



JiKJhk y„ JJri ^ z < " " - --inn,-- 

fc=i i=o 

Let (a n )„ e N be a sequence in R + and b > 0. Then, by standard results on the exponential rate of 
sums, we have 

n 

lim sup^ log a,kb n ~ k — max lim sup ^ log a„ , log 6 ^ . 

k — 1 
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Thus 



lim MogP Z0 (Z„ = z ) 

n~ too 

n 

< max { lim sup ± logE[Q„(z ) J[ fl(fi, n (0))] ; log E [Q(l) Za ] }. 



i=0 



We now prove that the first term always realizes the maximum. Using that /o(/o,n(0)) = z o/o°„ 1 (0) 
z Pi(Z„ = OlSY " 1 and fi{fi, n (0)) > f(0), we have 



E 



Qn(^0)ll^( ) 



i=0 



> zoE 



> zoE 



Q n {z a )¥ 1 {Z n = Q\£) Za - l J{Q i {l) 

Q„(^ )(<3n(o) n g 4 (i)) zo ~ n Qi(i) 



i=i 



> z o E[Q„(z o )g il (0)^- 1 jE[Q(l)^]", 
By definition of z , E[Q n {z )Q n (0) z °^ 1 } > 0, and we can conclude 

lim ilogP 20 (Z„ = zq) < limsupilogEtQ^zo)^^ - 1 ]^/^/,,™^))] 



i=i 



We end up the proof by checking the convergence of the sequence on the right-hand side above. 
We use again (3.4 1 and Lemma 3.1 to write 



b n : = P 20 + . . . + Z^ = Q,Z n = z Q ) 
= E 



E 



P Z0 (Z(°) + . . . + Zi n -V = \£)V(Y n =2b-l) 

n 

Qn(20)II/i(A»( )) 



(4.3) 



It is the probability of having zo-many individuals in generation n, where all individuals in gener- 
ation n have a common ancestor in generation n — 1. By Markov property, for /c = 1, . . . , n 

¥ Z0 (Z^ + + = 0,Z n = z o ) 

> ¥ Z0 (Z^ + ... + Zp-V = 0, Z k = z )F Zo (Z^l k +... + Z { ::*- X) = 0, Z n _ k = z ). 



The same subadditivity arguments as in the proof of Lemma |4. 1 1 applied to <p n yield the existence 
of the limit of — log <j) n and 

n 

lim ilogP Z0 (Z„ = z ) = lim ilog0 n = lim ± logE [o„(z )/ 0! „(0) Zo " 1 TT /• (/i,„(0))l (4.4) 

i=l 

This ends up the proof. 

□ 



5 Proof of Proposition |2.2 

5.1 Preliminaries on random walks 

In this section, we will shortly present some standard results on random walks which we use. In 
all following results, we take Sq = 0. We assume that there exists s > such that E[Xe" sX ] = 0. 
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This suggests to change to a measure P, defined by 

, , e- sx P{X G dx) 
P(X G dx) := 

where /i := E[e _sX ]. We note that ~E[X] = fi^ 1 E[Xe^ sX ] = 0, and that S is a recurrent random 
walk under P. In the following proofs, we use the change of measure described here. We define 

L n := min{5i,. . . , S n } 

and 

M n ;=min{Si,...,S n }. 

The following result is directly derived from [2] [Proposition 2.1]. For lattice random walks, the 
claims result e.g. from |34) [Theorem 6]. 

Lemma 5.1. Assume that E[X] = and V&r(X) < oo. Then for every 9 > 0, there exists 
d = d(0) such that 

3 

E[e~ es ";L n >0] ~ d n~2 (n -t oo) 

and 

E[e 9S ";M n < 0] - d n"5 (n -> oo). 
The following lemma results from [2] [Equation (2.5) therein]. 

Lemma 5.2. Assume that E[X] = and Var(X) < oo. TTien /or every c > Zarge enough, there 
exists d — d(c) such that 

3 

P(L„ > 0, S n <c)~dn2 (n -> oo). 

Remark: In the non-lattice case, the result holds for every c > 0. In the lattice case, c must be 
chosen such that P(0 < Si < c) > 0. 

From the previous results, it follows that 

Corollary 5.3. Assume that E[X] = and Var(X) < oo. Then for every 8 > 1, 

E[e- s " +eL »] =0(n- 3 / 2 ). 
Proof. We use a decomposition according to the first minimum of the random walk, i.e. let 

r„ := min{k G {0, . . . , n} | S^ = L n A 0}. 
Decomposing at the first minimum and using duality yields 

n 

E [ e -S„+*L„] = £ E [ e -(S„-fl5 fc ); T „ = fc] 

n 

= VE[e-^-^e( fl - 1 ) s *;7) b = A, min {S n - S,} > 0] 

^ — ' i—k,...,n 
k=0 



= STE[ e (0-i)S«. Tk = k ]E[e- s ^;L n _ k > 0] 

fe=0 

= £ E [e^ 9 - 1 ^ ■ M k < 0] E [e- s - fc ; L„_ fe > 0] . 



fc=0 
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The last step follows from duality (see e.g. [I]). Recall that by assumption, 6 — 1 > 0. Applying 
Lemma |5.1| now yields that there is a c < oo such that for n large enough 

E [ e S„+0L„] < c (_ + ^ / _ M3/2 fc 3/ 2 ) 



n 3/2 ( n _ fc )3/2 fc 3/2, 

3/2n 



A-l 

/ 2d 2d \ _ 

" C Vn372 + Ln/2j3/ 2 ^L«/ o 2J ^/J ~ 

which is the claim of the corollary. □ 

In the next lemma, we will use probability measures P + and P which e.g. have been introduced 
in [2J. Here, we recall the definition. Define the renewal functions u : M — > K and v : K — > K by 



= 1 + ^P(-S fc < x,max{5i,... ,5 fe } < 0), x > 0, 
fc=l 

oo 

u(a:) := l + 5^P(-S fc >aj,min{5i,... ) S , jfc }>0), x<0 

k=l 

v(0) = w(0) = 1, 
and elsewhere. Using the identities 

E[m(x + X);X + x>0] = x > 0, 

E[u(x + X);X + x < 0] = x < 0, 



(5.1) 



which hold for every oscillating random walk, one can construct probability measures P + and P _ : 
Define the filtration J- = (J r n ) n >o, where T n = o~(Q\, . . . , Q n , Zq, . . . , Z n ) . Then S is adapted 
to T and X n+ \ (as well as Q n +\) is independent of T n for all n > 0. Now, for every bounded, 
J-" n -measurable random variable R n we can define 

E+[i?„] = -^E x [R n u{S n );L n >0], x > 0, 

E-[i?„] = -^E x [R n v(S n );M n <0], x<0. 
v{x) 

The probability measures P+ and P~ correspond to conditioning the random walk S on not to 
enter (— oo,0) and [0, oo) respectively. 

More precisely, if R n —> R a.s. with respect to P + (resp. P~), then 

Km P(R n G -\L„ > 0) -> P+(i? G •) 

n— Voo 

lim P(i?„ G -|M„ < 0) -> P"(i? G ■)■ 

n— foo 

The first result is proved in [3] [Lemma 2.5] in the more general setting of random walks in the 
domain of attraction of a stable law. The proof of the second claim is analogous. 
We end up with an asymptotic result in the critical case, which is stated and proved only in the 
non-lattice case. 

Lemma 5.4. We assume that ELY] = 0, Var(Y) < oo, and that E[(log + £q(a)) 2+£ ] < oo for 
some e > 0. Then, for every c > 0, 

liminf P(Z„ > Q\L n > 0,S n < c) > 0. 
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Proof. The proof follows essentially J2j and we just present the main steps. First, 
Lemmas |5.1| and |5.2| ensure that for all 9,c > large enough, there exists d > such that 

E[e- es ";L n > 0] ~ d P(L n >0,S n < c) (n oo). 

Secondly, we recall the well-known estimate (see e.g. [Fj [Lemma 2]) 



(5.2) 



P(Z n > | £) > 



where % := E^li v{v ~ l)Q*(y)/ m Q 4 - Then, we rewrite 
1 



a.s., 



E 



e S " +E"=o lr ?i+ie s ' 
> E 



; L n >0,S n <c 
1 



> E 



; L n > 



= E 

> e- c/2 E 

,|n/2J 



(c-5 n )+ Al 

ip(U n ,V n ,S n );L n > 

fir " / V(C»,^n,5 n );L n >ol ) 



, -^n > 0, S n < C 



where U n := El=o ^+i e S *> ^™ := E™=L«/2j+i ^+i e5L " /2J Si and = (1 + u + u) 1 {c - 

e- c / 2 E 



z)+ A 1. Using (5.2), it becomes (with 9 = |) 



liminf P(Z„ > 0|L„ > 0, 5„ < c) > cT 1 liminf 



e- s "/ 2 <p{U n ,V n ,S n y,L n >0 
E[e-^/2;X„ >0] 



lim 



U n and 



Due to monotonicity and Lemma 3.1 in [2], the limits of Uoo 
Voc — limn^oo Ei=o^ exist and are finite respectively under the probabilities P + -a.s. and 

P _ -a.s. Thus all conditions of Proposition 2.5 in [5] are met. Applying this proposition with 
8 = 1/2, there exists a non zero measure v\/2 on K + which gives the convergence of the right-hand 
side above and 

liminfP(^ n > 0\L n > 0,S n < c) > / y(u,v, -z)P + (U 0o G du)P~(Voo G dv)v l/2 {dz) > 0. 

Note that in the function </?, z is changed to — z for duality reasons (see [5] for details). As 
and are a.s. finite with respect to the corresponding measures, this yields the claim. □ 

Remark. The proof may be adapted to the lattice case, by proving for example that 

n/2 

liminf P + (y^?7 i+1 e~ s * < d,S n /y/n G (a, b)) > 0. 

i=0 

We note that P + (Ei=o Vi+i e ~ Si < °°) = P~(Ei=i Vi eSi < °°) = 1 has been proved in [3] also for 
the non-lattice case. But the main remaining problem is that Proposition 2.5 in [2] is only stated 
for non-lattice random walks. The generalization of this result is a technically involved task and 
beyond the scope of this paper. 
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5.2 Proof of Proposition [272] (i) : p > 

Under Assumption [T] , we now prove that p > 0. It means that the probability of staying small 
but alive is exponentially small. The proof relies again on the Geiger construction and results of 
the previous section. 

We assume that there exists 7 > such that Q(0) < 1 — 7 a.s. and E[|A|] < 00. Let zq be the 
smallest element in X. Using (4.3) and ( |4.4[ ), we get 

lim i logP Z0 (Z„ = zo) = lim ± logE[P 2o (Z(°> + . . . + Z^ = 0, % = *> - 1|£ )P, (^ > 0|f )] 

n— ^00 n— foo 

n-1 

< liminf i log if IT P(£#) = 0|£) 

n-1 

= liminf i logfifexp f V logP(Z^ ) = 0|£) 

n— >oc n L V — ' 

3=0 

The fact that log(x) < x — 1 yields 



E 



7t — J. 

exp(^logP(Z« =0|£) 



j=0 



< E 



exp(-]TP(Z« >0|£) 



It remains to prove that this last expectation decreases exponentially. From (3.5 1, we get 

P(Z#> = 0|f) = -^/K/i,n(0)) = r^lr /j^W) = fci(/ iin (0)), (5.3) 
Pj — l,n 1 Jj\Jj,n\y)) 

where for s g [0, 1), 



/'(*) 



g(s) l-.s 
We will now show that <?(!) = /'(l), /i(l) = 1. As already noticed in [TJ], for every s £ [0, 1) 



"l _,*/<*> (0) 
5(S) = ^1- S fc! 



fc=0 



fc=i 



, fc -i/ fc) (0) 
J k\ ' 



Thus /'(l) = 3(1) = 1 and h(l) = 1. Moreover, /'(0) 7^ 1 ensures that for every k > 1 and s < 1, 
ks k ~ 1 < (1 + s + s 2 + . . . + s x ), so /i(s) < 1. A straightforward calculation shows that h has 
exactly one minimum in some Sg £ (0, 1). Adding that h is increasing for s > s and h(s ) < h(0), 
we have every t £ (0, 1) that 



h(s) < max {h(G),h(t)} for s < t, 



(5.4) 



First, we deal with fj. n (0) = V(Z n = 0\£, Zj = 1) in ( |5.3| . For this purpose, we use a truncation 
argument. Let a £ N be fixed for the moment and introduce 

QU) :=QV),l<j<a, Q(a) = Q([a, 00)). 

We refrain from indicating the dependence on a in our notation. The corresponding truncated 
random variables are denoted similarly, e.g. by X, S, /. Note that /"(l) < a 2 . By dominated 
convergence, we get that 

lim ELY] = ELY] > 0. 

a— J-oo 

Thus if a is chosen large enough, S is still a random walk with positive drift and E[/'(l)] > 1. Also 
note that with respect to the truncated offspring distributions, Z n is stochastically smaller than 
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Z n and thus P(Z n = 0) > F(Z n = 0). Applying this together with a well-known formula for the 
extinction probability (see e.g. [5] [Lemma 2]), we get that 

P(Z n = 0\£, Z 3 = 1) < Fi(Z n = 0\£) < 1 = — — ^ — — 

e -b n + J2 k=j rik+ie-^"-^' 

1 

- ~~ „2V°° „-2X fc + 1 -(S fc -S 3 ) ' 

where we used that fj = /"(1)//'(1) 2 < a 2 e~ 2X a.s. We now aim at bounding YlkLj ex P( — — 
(Sk — Sj)). First, the assumption Q(0) < 1 — 7 implies that X > log(7) a.s. Thus 

¥(Z n = 0\S, Z j = 1) < Px(Z n = 0\S) < 1 - ^ 

Next, we introduce the random walk S n := S n — en with < e < E[X]. It is still a random walk 
with positive drift and we have 

/in(0) < 1 1 = ■ 

Let us now consider the prospective minima (see e.g. [4] [p. 661]) of S which are defined by z/(0) := 
and 

v{j) := mf{n > v(j - 1) : S k > S n Vfc > n}. 
Then we can estimate for j > 1 (note that Sk > S v (j) f° r au & > Ki)) i > 1) 



a 2 ££=„y) e-(^- s -0'))e-( fe -"0-)) e 
<1 7 2 _ 1 7 2 (l-e' £ ) 



From (5.4), setting d := 7 (1 — e e )/a 2 e (0, 1), we get for j > 1, 

P(Z<"Ci)) > |f ) = 1 - /i Kj -)(/„( i)in (Q)) > 1 - max{/i(0),/i Kj) (l - d)} 

= min{l - h(0), 1 - h„y)(l - d)} =: 

From the classical random walk theory, Uj := f(j) — — 1) (and also Q u (j)) are i.i.d. random 
variables (see [2]). We now prove that for 5 > small enough that the probability that there 
are less than {<5n}-many prospective minima is exponentially small. Note that E[X] > implies 
E[i/(1)] < 00. Let < S < E[z/(1)]- 1 . Then 

P(tt{j > : u(j) < n} < Sn) < ¥{v([Sn\) > n) < p( U J > § n5 J ^ e~ 5n ^ d ^ , (5.5) 

j'=i 

where \& is the rate function of the process (J^ . —1 Uj) n , which is a random walk with nonnegative 
increments. Thus it remains to prove that ^(0) > for some 8 > E[f7i] = E[z/(1)]. From large 
deviations theory, we just need to check that the tail of v{l) decreases exponentially. This follows 
from 

00 

P(i/(1) > k) < P{Sj < for some j > k) < ^P(£j < 0) 

j=k 

00 „ p -A(0)fe 

< V e" A( °W = — — - 

-f-- 1 - e- A (°) ' 

j=k 
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where A is the rate function of S. This rate function is proper since log(l — 7) < X < a a.s. 
Adding that ELY] > ensures that A(0) > and we can conclude that *S?(0) > 0. 
Finally, we use that A v ij\ are independent to get 

n-l 

¥ Zo (Z n = zo) < E[exp ( - V ^n ] > \ S )) 

3=0 



< P(tt{j > : v(j) <n}< Sri) + E 



[5n\ 

exp i - J2 A »u) 

< P(tt{j > : v{j) <n}< Sn) +E[exp ( - A„ (1) )] ^ 



Recalling that A^) > and F(A v(j) > 0) = P(/i(0) < 1) = P(/'(0) ^ 1) > 0, we get 



E 



cxp ( - A v{j) 



< 1. 



Then (5.5) ensures that p > 0. 



□ 



Remark. To get Proposition 2.2 (i), Assumption^ can be replaced by assuming that there exists 
c such that r\ < c a.s. The proof is very similar. In this case, the truncation is not required and 
we may estimate 

1 1 



F(Z n = 0\£,Zj 



1 < 



0) < 1 



< 1 



5.3 Proof of Proposition 2.2 (ii) : g < A(0) 



Here, we prove the second part of Proposition |2.2| which ensures that g < A(0). It means that 
small but positive values can always be realized by a suitable exceptional environment, which is 



'critical'. We focus on the nontrivial case when A(0) < 00. The proof of Proposition 2.2 (ii) can 



then be splitted into two subcases, which correspond to the two following propositions. 
Proposition 5.5. Under Assumption^ and P(Y < 0) > 0, we have p < A(0). 
Proposition 5.6. Assume that P(Y > 0) = 1 and P(Y = 0) > 0. Then 

p< -logP(Z = Q) =A(0). (5.6) 



Proof o] Proposition^^ We recall that X := {j > 1 : F(Q(j) > 0,Q(0) > 0) > 0}. We use a 

standard approximation argument and consider the event E x<n :— {min; = i „ Xi > x} for x < 0. 

Then, P(Y > x) > since we are in the supercritical regime and for every s > 0, E[|A|e~ sX |A" > 
x] < 00. As P(Y < 0) > 0, we may choose x small enough such that P(a; < X < 0) > 0. Then 
E[|Y|e _sX |Y > x] tends to infinity as s — > 00. Moreover E[e~ sX |Y > x] is differentiable with 
respect to s for s > 0. We call s — v x a point where the minimum of this function is reached. In 
particular, 

inf K[e~ sX \X > x] = E[er v * x \X > -x], £E[ e - sX \X > x]\ s=v = E[Xe-^ x \X > x] = 0. 



The second part of Lemma 4.1 ensures that for every z G X and for every sequence k n = o(n), 
lim A logP 2 (Z„ =z)= lim I IogP,(l < Z n < k n ) = -p. 
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Let us now change to the measure P, defined by 

e -»*vf(X G dy\X > x) 



P(* G dy) 



(5.7) 



where /x := E[e \X > x\. Under P, E[X] = and S is a recurrent random walk. 
Let c > be so large such that P(L„ > 0, S n < c) > for every n. Then 



,(1 < Z n < k n \£);L n > 0,S„ < c 



(5.8) 



We note that P z (l < Z n < k n \£) = P z (Z n > 0\£) - P z (Z n > k n \£) a.s. and by Markov inequality, 



,{Z n > k\£) < 



a.s. It ensures that 



s (l < Z n < k n \E x , n ) > M "E P(Z„ > 0\£) - ze c /k n ;L n > 0,S n < c 



Plugging this into (5.8) and setting b n :— P(L„ > 0, S n < c), we get 

P,(l < Z n < k n ) >¥ z {l<Z n < k n \E Xin )¥(E x<n ) 

= n n b n [p(Z n > 0\L n > 0, S n < c) - ze c /k n ] P(X > x) n . (5.9) 

By construction of E, Var(X) < fj,~ 1 E[X 2 e~ VxX \X > x] < oo. Then from Lemma \5/2\ we have 
b n = 0(n~ 3 / 2 ) and limn^oo ilog&„ = 0. Let k n = rT 1 / 2 . The fact that Assumption [2] holds 
under P entails that it holds also under P. Indeed, 

E[(log+£ Q (a)) 2+£ |X > x] = M E[ e ^(log+£ Q (a)) 2+£ ] 

>/i e ^E[(log + £ Q (a)) 2+e ], 



as X > x (x < 0) P-a.s. Thus we can use Lemma 5.4 and ( |5.9[ ) to get 

liminf i log'P 2 (l <Z n < k n ) > log/x + logP(A > x) 

= logE[e~^ x |A: > x] + log¥(X > x) 
= -sup{ - logE[e" sX ;A: > x]}. 

s<0 

By monotone convergence, we let x — > -co and 

liminf i log P,(l < Z„ < fc„) > -sup{ -logE[e- sX ]} = -A(0). 
n-K» n s < 

As k n = o(n), we apply Lemma |4.1| to end up the proof. 



□ 



Proof of Proposition^^ As ¥{X > 0) = 1, we have P(S„ = 0) = ¥(X = 0)" and A(0) = 
-logP(X = 0). 

If P(Q(1) = l\X = 0) = 1, the proof is trivial. So let us work with the assumptions P(X = 0) > 
and P(Q(1) = 1\X = 0) < 1. 

By conditioning on the environment, we get for z£l that 



P z {Z n =z)> ¥{X = 0)" • ¥ z (Z n = z\X x 



0, 



0). 



For simplicity, we introduce a new measure P on the space of all probability measures on No with 
expectation 1. It is defined for every measurable A C A by 

G A- vriQ = 1) P(Q G A] vtiq = 1) 



P(Q G A) 



P(m Q = 1) 



¥(X = 0) 
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4.1 



Note that P(X = 0) = 1 and there exists z > 1 such that P(Q(z) > 0, Q(0) > 0) > 0. With 
respect to P, (Z n : n £ No) is still a branching process in random environment. By Lemma 
there exists p £ [0, oo) such that 

p = - hm i \ogP z (Z n = z) 



lim i logP,(Z„ = = 0, . . . , X n = 0) 

n— >oo 

n-l 

lim IlogEb^z)/*- 1 ^) J] //(/^(O)) 



= 0, . . . , X n = 



Next, we use convexity arguments. First, for all i < k and s £ [0, 1], /j,fe(s) > 1 — — s). 

As P(/£ fc (l) = 1) = 1, we get 



/i,fc(s) > s P-a.s. 



(5.10) 



Also recall that / , n (0) = / 0j „-i(/„(0)) and by ( pO] ), /,,„(0) > / n (0) = Q n (0). Thus, for every 
a £ N, 



n-l 



hminfilogE Q n (z)ft-\0) ]J tf(/<,„(°)) 



A! = 0,...,X„ = 



n-l 



>liminfilogE Qn^QnCO)*" 1 IJ /i(/n-»,«(0)) II /<(A™(°)) 



z=n— a+1 



For e > fixed, we choose fc = k £ £ N large enough such that P(Q([1, k]) > e) > 1 — e. Then, 
conditionally on {Q([l,fc]) > s}, f'(s) > J2j=iQ( k ) sk > £sk a - s - for s G [0,1]. Using this 
inequality, (5.10) and /i, n (0) < /i in _i(0), we have 

n — a n— 1 



E 



> E 



Qn(z)Qn(0) Z - 1 I] £(/«-«,n(0)) I] Gl, ' " ' > Qn-a 

i— 1 i— n— a+1 

n—a 

Q„(^)Qn(0) 2 - 1 n#(/"-«.n(0)) 

2 = 1 

n-l 

x 1] fi{Qn(0));Qn-a+i([hk])>e,...,Q n - 1 ([l,k})>e 

i—n—a-\-l 

n—a 

QnWQniO)'- 1 n/i(/™-l(0)) 



(^i , . . . , Qn—a 



> E 



;=i 



x J] eQn(0) fc ;Q„_„+i([l,fc]) >£,..., Q„_i([l,fc]) >e 

z— n— a+1 

n—a 

> E[ J] /;(/„- a ,„-i(0)); Q„_ 0+ i([l, fc]) >e,... , Q„-i([l, fc]) > e 



2=1 



1 ) * * • ) Qn—a 
!)*••) Qn—a 



x E 



s a - 2 Q n (z)Q n (0) z ~ 1Ha ~ 2)k 



where the second expectation is strictly positive as P(Q(z) > 0, Q(0) > 0) > 0. The product of two 
generating functions (and thus the product of finitely many) is again convex. Indeed generating 
functions, as well as all their derivatives are convex, nonnegative and nondecreasing functions, thus 



(f9)" = f"9 + 2g'f + fg">0. 
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Similarly, the product of the derivatives of generating functions is again convex. For more details 
on the product of nonnegative, convex and nondecreasing functions, we refer to |30) . Applying 
Jensen's inequality to the convex function U^Z^fl, the independence of the environments ensures 
that 



E 



Ql i • • • i Qn—a 



'{ II /i(/n-a,n-l(0))5 0«-a+l([l. k}) >£,... , Q n ^{[l, k]) > E 

1=1 

n—a 

> l[ #(E[/„_ a , n _i(0); Q„_ a+1 ([l, k}) > e, . . . , Q„_i([l, k]) > e\Q u ■ ■ ■ , Qn-a\) 

n—a 

= W /i( E [/o,a-i(0); Qi([l,fc]) >e,..., Q a _i([l,fc]) > e]) P-a.s. 

i=l 

Using this inequality yields 

n 

hminf i logEk'-^O) J] £ (/<,»(<>)) 



X! =0,...,X n = 
> hminf I log ( E[ J] //(E[/ Oia _i(0); Qi([l, fc]) >£,..., Q„_i([l, k}) > e\) 

'e a ' 2 Q n (z)Q n (0y~ 1 (0)Q n (0)^ k 



i=l 
n—a 



i=l 



X E 



= lim inf - log E 

n— >oc n 



logE 



/ / (E[/o,„-i(0);Qi([l,fc])>e,...,Q _i([l,*])>e]) 
/'(E[/ , o -i(0);Q 1 ([l,fc])>e,... ) Q a _ 1 ([l,fc])> £ ]) 



Finally, Z is a critical branching process in random environment under the probability P so 
P(Z a _i = 0|£) = /o,a-i(0) — > 1 P-a.s. as a — > oo (see e.g. 133'). Letting a — > oo, e — > 
and recalling that P(Q([1, k]) > e) > 1 — e yields by dominated convergence 



logE /'(E[/o, tt -i(0);Qi([l,fe])>e,...,Q«-i([l s fe])>e]) -> logE[/'(l)] = 0. 



Then, 



n-l 



hminf i logE /'"HO) [J /i(A»(0)) 



x 1 = o,...,x n = o 



> 0. 



This yields the claim. 



□ 



Remark. Note that the bound f'(s) < /'(l) for s £ [0,1] immediately yields that 
lim sup^^^ ^logP z (Z n = z) < logEf-X"]. In particular, we recover that for a BPRE with X = 
a.s., the probability of staying bounded but positive is not exponentially small. 



6 The linear fractional case : Proof of Corollary 2.3 



In this section, we assume that the offspring distributions have generating functions of linear 
fractional form, i.e. 

1 - s 



1 



m 1 +bm 2 (1 — s)/2' 
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where m = f(l) and b = /"(l). 

Under this assumption, direct calculations with generating functions are feasible, i.e. we can 
explicitly calculate the generating function of Z n , conditioned on the environment. We also assume 
that E[|X|] < oo, P(Zi = 0) > and that either F(X > 0) = 1 or Assumption [2] hold, such that 
Proposition 2.2 (ii) holds. 

In the next subsection, we prove Corollary |2.3| It gives an expression of g which depends on the 



sign of E[Xexp(— X)]. Afterwards, we prove Corollary 2.4 which concerns the MRCA. Let us 
dehne r]k := fe^m^ 2 /2 and recall that fj >n = o . . . o f n . Then for all n € N and s € [0, 1] (see 
El p. 156]) 

(1 - s) 



/o,n(s) — 1 _ , • 

e 5 - + (1 - s) }_^ k=0 r] k+1 e b >> 



resp. 



/j,n(0) 1 _fo _q.\ i v^™ _("<?, ,_<?.'!' (6-1) 

Let us state some direct consequences resulting from this formula which will be used later. Taking 
the derivative, 



/o nOO = i a- ( 6 - 2 ) 

' (e-^ + (l- S )E^- ^ 

Note that for every s£ [0, 1), applying (|6.1| 



C 1 - s )ELo 1? ?fc+ie Sfc 



/o. n W(l-,s) 2 = 



n2 e 



< e 



1 



e 



~ 5 " (1 - /o,«(0)) 2 = e- s "P(Z n > 0\£) 2 . (6.3) 



Moreover, 

f'.(a\ = 

{e-^ +71,(1-8)) 



and we can now compute the value of g. 



6.1 Determination of the value of g 



X] 



> 



By Proposition 2.2 (ii), p < A(0). Then it remains to prove that p = — logE[e _x ] if E[Ae 
and p > A(0) otherwise. For that purpose, we use the representation of p in terms of generating 
functions. Combining ( |6.1[ ) and (6.4 1 we get 



/j(/i,»(0))=e-*'( 



—X4 l — Xj 



-(Sn-Si) | v^" -(Si. , -S,0 



)" 



-X; / Z ^ k =3- 



-{Sk — l — Sj) 



p (z„>o|z J _ 1 = i,g; 

P(Z„ >0|Z 3 - = 
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Since P(Z n > 0\Z n = 1) = 1, we get 

n 

Pi(z„ = i) = Ex [JJ /<(/,•,„(())) 



n 

3=1 



_ x , P(Z n >0|Z J -_ 1 = l,g) 2 
P(Z„ > 0|Zj = 1,£) 2 



= e 



P(Z„ > 0|£,Z = 1)< 



(6.5) 



First, we consider the case E[Xe ] > 0. Bounding the probability above by 1 immediately yields 



n 

E [lI/j(/7>( ))] <E[e- 5 "]=E[e 



3=1 



„-Xm 



Thus p > — logE[e ]. To get the converse inequality, we change to the measure P, defined by 



P(X £ dx) 



e- x V{X £ dx) 



E[e 



-XI 



This measure is well-defined as E[A 2 e x ~\ < oo implies E[e x ] < oo. Then by Jensen's inequality, 



Ei 



P(Z„ > 0|£) 2 j = E[e- x ]"E[Pi(Z„ > 0\£) 2 ] > E[e- x ] n P 1 (Z n > 0) 2 . 



We observe that E[X] = E{Xe~ x ] > 0, such that under P, S n is a random walk with nonnegative 
drift. It ensures that the branching process is still critical or supercritical with respect to P. Thus, 
under Assumption^ and as E[X] = E[X 2 e~ x ] < oo, V(Z n > 0) > Cn"2 for some C > as n -> oo 
(see e.g. [1] for the critical case, whereas P(Z n > 0) has a positive limit in the supercritical case). 
Letting n — > oo and adding that 1 £ I, we get 

n 

p = - Urn i log E [ JT f. „ (0))] < - log E[e- X ] . 

Secondly, we consider E[Xe~ x ] < 0. Then there exists v £ (0, 1) such that ELYe -1 ^] = and 
we change to the measure P defined in ( |5.7| ) (here without truncation). Applying this change of 
measure and the well-known estimate P(Z n > 0\£) < e L " A0 a.s., we get that 



E 



(Z„ > 0\£, Z = l) 2 < E[e~ vX ] n E 



J-l+»)S n +2L n A0 



Note that L n A < min(5„, 0) and v £ (0, 1) that (-1 + v)S n + 2L„ A < 0, and thus 

„ >0\£,Z = 1) 2 } <E[e- x ] n . 



E 

This yields p > -logE[e- yX ] = A(0) since A(0) = sup s<0 {- logE[e sA ']} and the condit 
ElATe -1 ^] = implies that the supremum is taken in s 



ion 



6.2 Proof of the limit theorems for the MRCA 



The three cases (i-ii-iii) in Corollary 2.4 result respectively from Lemma 6.1 Lemmas 6.2 6.4 6.5 



and Lemma |6.7| below. In all this Section, we assume that the assumptions of Corollary |2.3| are 
met, i.e. that ELY 2 e~ x ] < oo, E[\X\] < oo and either ¥(X > 0) = 1 or Assumption [2] holds. 
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Lemma 6.1. IfE[Xe- x ] < 0, th 



liminf Pi (MRCAn = n\Z n = 2) e (0, 1) 



and 



liminfPi (MRCA n = l\Z n = 2) e (0,1). 

n— >oo 

Proof. Conditionally on £, the branching property holds and ensures that 

F 1 (MRCA n = n\Z n = 2) > P l{Zl = 2) . ^f^Z^ , 

which corresponds to two subtrees being founded by Z = 1 and staying equal to 1 in all generations. 
Next, we recall that linear fractional offspring distributions are stable with respect to compositions 
and have geometrically decaying probability weights (see e.g. [25]). In particular, 



and therefore 



\(Z n =2\£)<¥ 1 (Z n = 1\£) 

E[Pi(Z^i = l\£) 2 ] 



\{MRCA n = n\Z n = 2) > Pi(Z a = 2) 



i(Z n 



(6.6) 



(6.7) 



Note that E{Xe~ x ] < implies that there exists v € (0, 1) such that E\XeT vX \ = 0. Thus we 



can apply the same change of measure as in the proof of Proposition |5.5| With the definition of P 
therein (again without truncation), we get 

EpPi(Z„ = l\£) 2 } _ E[e" s »Pi(Z„ = 1\£) 2 } fe , 



\(Z„ = 1) 



From (6.5), we know that 

Pi(Z B = 1\£) = &J1) 



E[e" s »Pi(Z n = 1\£)] 



e - s "Pi(Z„ > 0\£) 2 < e 



h2L„ 



a.s. 



(6.9) 



Combining this with Jensen inequality yields for every c > 0, 

E^-PxCZ, = l\£) 2 } > E[e- aff » + " s »Pi(Z n > 0\£) 4 ;L n > 0,S n < c] 

> e^+^EfPi^ > 0\£) 4 :L n > 0,S n < c] 

> e^ 2+ ^ c P(L n > 0,S n < c)E[P x (Z n > 0\£) 4 \L n > 0,S n < c] 

> e ( - 2+1/)c P(L„ > 0,S n < c)P{Z n > 0\L n > 0,S n < c) 4 
>dP(L n >0,S n <c), (6.10) 



for some constant d > 0, where the last line follows from Lemma 5.4 For the denominator in (|6.8|) , 



by (6.9), we get similarly 

E[e" s »Pi(Z„ = 1\£)} < E[ e -(i-»)Sn+2L n ^ 

Finally, we use ELY 2 exp(— X)) to ensure that ELY 2 ] < oo and apply Lemmas 

\{Z, 



5.1 



and 



5.2 



Then, 



Efe 



vS n 



v • l|g) 2 ] ^ r . . P(L n >0,S n <c) ^ 
limmt —7 — „ — -; — — > limmf ^—7 — 77 , „ , „. . > 0. 

n->oo Efe^P^Z^ = 1|£ )] ~ n->oo E[ e -{.l-v)S n +2L n ] 



Recalling (6.7) yields the first part of the lemma. 
Similarly, we note that by the Markov property, 

P^MRCAn = 1, Z n = 2) > Pi(Z n _i = l)P(Zi = 2). 



Recalling that Pi(Z n _i = 1) > Pi(Z„_i = 2) from (6.6) yields the second claim. 



□ 
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The next three lemmas cover the 'intermediate' regime. First, we prove that the probability of 
{MRCA n — x n } doesn't decay exponentially for every 1 < x n < n. 

Lemma 6.2. IfE[Xe~ x } =0, then for every 1 < x n < n 

lim i logPi {MRCA n = x n \Z n = 2) = 0. 

n— >oo 

Proof. Let E[Xe~ x ] = and note that MRCA n > 1. Conditionally on £, the branching property 
holds and guarantees that as in the previous proof that for every 1 < x n < n, 

EpPi^ = 1\S ) 2 1 
¥ x (MRCA n = x n \Z n = 2) > P x (Z n _ SB = 2) • = 2) • 

As E[ATexp(— X)] > 0, we know from the previous subsection that logE[e _J<: ] = limn-j.oo Pi(Z n = 
2). Thus we get for 1 < x n < n, 

liminf ± log Pi (MRCA n = x n \Z n = 2) 

71— i-OC 

> hminf ((1 - x n /n) logE[ e - x ] + ± logE[Pi(Z aB = 1|£) 2 ] \ - logE[ e - x ]. 

For the last term, we use the change of measure of the previous lemma with v = 1, i.e. 

e~ a: P(X G da;) 



P(X G dx) 



E[e 



Applying (6.10), we get that 

E[Pi(Z In = 1|£) 2 ] > d E[e- x ] a "P(i x „ > 0, S Xn < c) 



As E[Ae ] = 0, S is a recurrent random walk under P. Using again Lemma 5.2 ~P(L n >0,S n < 
3 1 ' 

c) ~ (in 2 ; and 

liminf i log Pi (Afi?CM„ - x n \Z n = 2) 

71— >O0 

> logE[e~ x ] - logE[e~ x ] + liminf j - x„/n logE[e~ x ] + x„/nlogE[e" x ] 1=0. 



It gives the expected lower bound, whereas the upper bound follows is simply due to the fact that 
the probabilities are less than 1. □ 

The following lemma is required to prove limit results for the MRCA in the intermediately super- 
critical case. To avoid technicalities, we impose an additional moment condition. 

Lemma 6.3. Let E[Ae~ x ] = and 7 = E[e~ x ]. Assume that E[/"(l)e~ x ] < 00. Then 
< liminf 7 -"n 1/2 P X (Z„ = 2) < limsup7"n 1/2 Pi(Z n = 2) < 00, 

i.e. Vi{Z n = 2) is of the order ^n' 1 / 2 . 



Proof. Using the change of measure to P and (6.9), 



\{Z n = 2) =E[e- s -V 1 (Z n >0\£) 2 ] 

= 7 "E[P(Z„ > 0|£) 2 ] > j n P(Z n > 0). (6.11) 
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Under P, E[X] = and Z is a critical branching process in random environment. Under our 
assumptions, there is a constant d > such that (see [3] [Theorem 1.1]) 

P(Z n > 0) - d n~\. (6.12) 



Following the proof of the upper bound, using (6.11) and that P(Z n > 0) < e Ln a.s., we have 

Pi(Z„ = 2) = 7 n E[P(Z„ > 0|£) 2 ] < 7 "E[e 2i "]. 

Our assumptions imply E[X 2 ] < oo and thus, as a direct consequence of 0] [Lemma 2.1 and 
I-oo e x u(~x)dx < oo], 

E[e 2L "] = 0(n-^ 2 ). 

This proves the upper bound. □ 
The next lemma describes the probability of {M RC A n — n} in the intermediate regime. 
Lemma 6.4. Let E[Xe~ x ] = and assume thatE[f"(l)] < oo. Then 

< liminf n Pi (M RCA n = n\Z n = 2) < limsupn ¥ x (MRCA n = n\Z n = 2) < oo, 
i.e. ¥x{MRCA n = n\Z n = 2) is of the order rT 1 . 

Proof. First, the event {MRCA„ = n} implies that there are at least two individuals in generation 
1 and that from this generation on, at least two subtrees survive until generation n. We use the 
branching property and a decomposition according to the two subtrees which survive and get that 

P^MRCAn = n, Z n = 2\£) 

= £ Q) Pl ( Zl = ^) p ( z » = l \ £ > = l ? v (Zn = o\e, Zi = i) k - 2 
fc(fc-i) T 

< P(Z n = 1\£, Z, = l) 2 /r(l) a.s., (6.13) 



< ]T 2 F ^ = k \ £ )mn = l\£, Zt = 1) 

k=2 



since /"(l) = J2k?=2 k ( k ~ l)Pi(^i = k\£). Using (6.9) and independence yields 

¥ 1 (MRCA n =n,Z n = 2) 
< E[/"(l)]E[e- 2S - 1 P(Z„_ 1 > 0|£) 4 ] < E[/"(l)]E[e- 2S '' l - 1+4Z '"- 1 ]. (6.14) 

Again, we change to the measure P. Note that the assumptions of the lemma E[/"(l)] < oo 



ensures that Var(X) < oo. Using also Corollary 5.3 s we get 



Pi(MRCA n = n,Z n =2) < E[/"(l)] 7 "- 1 E[ e - s "- 1+4i "- 1 ] = 0(y i n~ 3/2 ). 
Inserting this and applying Lemma |6.3| yields 

limsupn Vi(MRCA n = n\Z n = 2) < oo. 

n— ►oo 

For the lower bound, we use similar arguments. First, 

Pi (MRCA n = n, Z n = 2\£) > V 1 (Z 1 = 2\£ )Pi(Z n -i = l\£, Z x = l) 2 a.s. 
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Let c > 0. Taking the expectation and using (6.9) yields 

Pi (MRCA n = n,Z n = 2) 

> ¥ 1 {Z 1 = 2)E[e- 2S - 1 P 1 (Z n _ 1 > 0\£) 4 } 

> P x (Zi - 2) 7 "- 1 e- c E[P 1 (Z„_ 1 > 0|£) 4 |£ n _i > Q,S n -i < c]P(L n _i > 0,S n -i < c). 



Moreover, by Lemma 5.4 and Jensen's inequality, 

liminf E[P X (Z„ > 0|£) 4 |L„ > 0,S„ < c] > 0. 

n— foo 

Applying Lemma |6 . 3| again , we get that 

liminf n Pi (MRCA n = n\Z n = 2) > 0. 

n.— >co 

□ 

The next Lemma describes the probability that the MRCA is neither at the beginning nor at the 
end: 

Lemma 6.5. Let E[Xe~ x ] = and E[/"(l)/(l - /(0)) 2 ] < oo. Then for every S G (0, 1) 
< liminf n 3/2 P x (M RCA n = \Sn]\Z n = 2) < limsupn 3/2 F^MRCAn = \8n\\Z n = 2) < oo, 

n >oo n — ^.qq 

i.e. Pi {MRCAn = \5n\\Z n = 2) is of the order n~ 3 / 2 . 

Proof. First, the event {MRCA n = \Sn\,Z n = 2} implies that the two individuals in generation 
n stem from one individual in generation n — \Sn~\ = |_(1 — 5)n\. If there are k individuals in this 
generation, there are k possibilities for the ancestor from which the two surviving individuals in 
generation n stem from. All others have to become extinct. We use the branching property and a 
decomposition according to the two subtrees which survive to get that a.s. 

¥{MRCA n = \5n\,Z n = 2\£) 

OO 

= X>(%i-*)»J - k\£)kP(Z n = 0\£, Z l(1 _ s)nj = h- 1) 
fc=i 

V(MRCA n = \5n],Z n = 2\£ , %i_ 5)nJ = 1) 

CO 

= 5Z P (%i-<5)«J = k\£)k¥{Z n = 0\£,Z l{1 _ s)ni = l)* 5 " 1 

k=l 

F(MRCA n = {5n],Z n = 2\£, Z L(1 _ 4)nJ = 1) 



Next, we set 



and note that 



Thus we get that a.s. 



8 := (s(£) = F(Z n = 0\£, Z L(1 _ ff)nJ+1 = 1) 

P(^n = 0|£ , ^L(l-«5)nJ = 1) = /L(l-<5)nJ ( S )- 



P(Mi?C^„ = r<5n]^« = 2|5) 
= fua-S)n}(f[(i-s)n}(s))V 1 (MRCA n = \6n],Z n = 2\£,Z [{1 _ s)ni = 1) (6.15) 
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Next, using (6.13) and (6.9), we get that 



\{MRCA n = \Sn] , Z n = 2\S, Z l{1 _ s)nl = 1) (6.16) 

< /ftl-i)nj+l(l) P (^» = l\£,Z l(1 - S)nl+1 = l) 2 

= /l(i- A > iJ+ i(1)^ 2(5 ^ Sl(1 - 5) " j+i) P(^ > 0\S,Z [(1 _ s)ni+1 = l) 4 

= /r (1 - 5H+ i(l)e- 2(S "- SL(I - 5 >"J +l) (l ~ «) 4 . (6.17) 



Combining (6.15) and (6.17), we have a.s. 

¥(MRCA n = \5n],Z n = 2\£) = /^ L(1 _ 5) „j (/ L( i- 5) „j ( S ))/[' (1 _ 5)nJ+1 (l)e- 2(S "- s L( 1 -)"J +1 



Hi-*) 4 



Moreover, from (6.3) 



/o.L(i-^j(/L(i-5)«j(s))(l-/L(i-^j(^)) 2 < e- s ^-^iP(Z l(1 _ s)ni > 0\£) 2 a.s. 



So 



¥{MRCA n = [8n] , Z n = 2\£) 



'L(i-«)"J 



(!-•)' 



As already used in [12| [Proof of Lemma 1] , we have for a generating function / of a random variable 
R 



l-s 



^Vp(i? > fc), 



fc=0 



which is obviously increasing in s. Thus we get 
1 - s 



< {nZ V(1 - S )n\ > 0|Z L( i_tf)„j_i = !,£))' 



l-/ L (i-*)*jOO " -^-o,^ i-/ L(1 _ 5)nJ ( )- 

Combining these identities and using the independence of the environments yields : 

P(MRCA n = \Sn\ , Z n = 2) 
<E[e- s L(i-*)»jP(Z L(1 _ 5)nJ >0|f) 2 ]E[/"(l)/(l-/(0)) 2 ]E[e- 2S r-i-P 1 (%„ 1 _ 1 >0|f) 2 ], 

where we recall that by assumption E[/"(l)/(l - /(0)) 2 ] < oo. As we have proved before, for 
every S £ (0, 1), 

E[e- fi, " 1 -*)»JP(Z L( i_ 3)nJ >0|£) 2 ] < 7 L(i-- 5 )«JE[e 2L L<i-^)"j] = 7 L(i-«)"Jo(n -1/a ) 

and 

Ele-^-i-ip^Z^.! >0|£) 2 ] < 7 ^l-i E [e- s ^i- 1+2i ^i- 1 ] = ^"^(n" 3 / 2 ). 
Together with Lemma |6.3[ this yields the expected upper bound: 

limsupn 3/2 Pi (MRCA n = \Sn] \Z n = 2) < oo. 

n— foo 

For the lower bound, we use 

Pi (M#CA„ = [Sn\,Z n = 2\£) 
> Pi(Z L( i_ tf)nJ - l|f)Pi (Mi?CA„ = r<5nl,Z n = 2|f , Z L(1 _ 5) „j = 1) a.s. 
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Both terms are independent and from Lemma 6.4 we get 

liminf 7 -L' 5 «J+ 1 n 3/2 Pi ( Mi?c , An = r §n i Zn = 2\Z Ul _ s)n] = 1) > 0. 

n— >oo 

From the previous lemmas, 

liminf 7 -L(i-^™J n i/2 Pi(z = 1) > 0. 



Thanks to Lemma 6.3 we obtain the expected lower bound. 
For the next proof, we require the following auxiliary result. 
Lemma 6.6. We assume that E[Xe~ x ] > 0. Then for every c > 0, 



□ 



lim I log inf F z (Z n = 2) = -oo, 

n— too z>cn 2 



lim sup i log Pi (1 < Z n < cn 2 ) = logE[e" x ]. 



Proof. First, we observe that on the event {Z n = 2}, at most two initial subtrees survive until 
generation n. Using the the branching property conditionally on £, we have a.s. 

inf ¥ z (Z n = 2\£)< V ( )Pi{Z n = 0|£) Lc " 2j ~ 2 Pi(l < Z n < 2\£) 2 

— i — — 2 V / 



k=[cn 2 ] 



+ J2 ( i m(Zn = 0\£)t cn ^ 1 ¥ 1 (l<Z n <2\£). 

k=[cn 2 i ^ ' 

Again, we use the geometric form of LF distributions to see that P(l < Z n < 2\£) < 2P(Z n = 1\£) 
a.s. Next, we use 



^k{k-\)a k - 2 < 



a 



ri-2 00 



{I -a) 3 ' 



£ k a k - x < 



to get 



inf F z {Z n =2\£) 

z>cn 2 



< n 1 



h(Z n = 0\£) 



Lcn 2 J-2 ]] 



(l-Pi(Z„ = 0|£)) 3 
MZ„ = 0|£)L™ 2 J- 2 



Y {Z n = \\£Y Pi(Z n =0|f) 
V c — 



[cn 2 ]-lT, 



l(Z n = l\£) 



1 {Z n = l\£) 2 P 1 (Z n = 0\£) 

h C 



(l-Pi(Z„ = 0|£)) 2 



Pi(Z„>0|£) 3 Pi(Z„>0|£) 2 
c 2 Pi(Z„ = 0|£) Lc,l2j - 2 e- 2S "Pi(Z„ > 0|£) + cPi(Z n = 0|f)L cna J- 1 e -' s " 



by using (6.5). Finally, as P(Z n > 0|£) < e 5 ™ a.s., we get that a.s. 



inf F z (Z n = 2\£) < (c 2 + c)n 2 Y x {Z n = 0\£) 

z>cn 2 



Next, we use again the change of measure 



P(le dx) 



e~*F(X € dx) 



E[e 



-X] 



Then 



inf, V z (Z n = 2) < (c 2 + c)n 2 E[e- x ]"E[P(Z co = 0|£)L C " J" 2 ]. 

z>cn 2 
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Using Jensen's inequality yields 

limsup i log inf V z {Z n = 2) < logE[e~ x ] + limsupnE[logP(Z 00 = 0\S)]. 

n-^-oo z>cn 2 

Finally, note that E[X e - x ] > implies ELY] > 0. Thus P(P(Z oc = 0|£) < 1) > and therefore 
E[logP(Z 00 = 0\£)] < 0. This yields the first result. 

For the second claim, we use again the geometric form of the probabilities of of LF distributions 
to get 

Pi(l < Z n < cn 2 \£) < \cn 2l [ ¥ 1 {Z n = l\£). 
Thus, taking expectations yields 

limsup A logP^l < Z n < cn 2 ) < limsup A log (Lc^P^ = 1)) 

= limsup i logPi(Z„ = 1) = logE[e- x ] . 
where the last result has been shown in the previous subsection. □ 
Lemma 6.7. We assume that E[Xe~ x ] > 0. Then, for every S G (0, 1], 

limsup i logF^MRCAn > Sn\Z„ = 2) < 0. 

n— ¥ oo 

Proof. First, we recall that the event {Mi?CA„ = [Sn\ +1} implies that there are at least two 
individuals in generation n — [6n\ and that from this generation on, at least two subtrees survive 
until generation n. As in the preceding lemmas, we use the branching property and a decomposition 
according to the two subtrees which survive and get that 

Pi (MRCA n = I5n\ + 1, Z n = 2) 

< J2 ) p i( z n-[8 ni = fc)E[P 1 (Z L5 „ J = l\£) 2 ¥(Z [Sni = 0\£) k - 2 ] + inf ¥ z (Z [Sni = 2) 
<A n + B n , 

where 

A n := n 4 Px(l < Z n _ [Sni < n 2 )E [P : {Z [Sni = 1\£) 2 ], B n := inf P z (%„j = 2). 
Letting n go to oo and applying Lemma |6. 6| yields 

limsup i \og¥x{MRCA n = [5n\ +l,Z n = 2) 

n— >oo 

< max < lim sup ~ log A n , lim sup ~ log B n > 

^ n— *-oo n— >oc ' 

= max < limsup ^ log A n , — oo ^ = limsup ^ log A n . (6.18) 
Next, let us treat the term named A n . By Lemma |6.6| 

limsup ilogP(l < Z n _ [5n] < n 2 ) = (1 - 5) logE[e- x ]. (6.19) 



Using ([6J5} and ¥(Z n > 0\£) < e L " a.s. yields 

E[Pi(Z LanJ = l\£ ) 2 ] < Efc-^-jp^Z^j > 0|£) 4 ] < E[e- 25 L^j+"L^j]. ( 6 .20) 
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We recall the change of measure 

Then E[Xe~ x ] > assures that E[X] > and under P, S is a random walk with positive drift. 



From \6.18\ , ( |6.19[ ) and ( |6.20| ) we get 

limsup i log Pi (M#CA„ = [5nJ + 1, Z n = 2) 

n—^oo 

< (1 - <5)logE[e~ x ] + limsup i log E[e _2S Li"J+ 4i H"J] 

n— foo 

= (1-5) log E[e~ x ] +<5 log E[e~ x ] + limsup A log E[e _ ' s l*«J+ 4Z 'L««J ] 
= log E[e~ x ] + lim sup A log E[e" s L^J + 4L Ls»j ] . 

As E[Ae~ x ] > 0, we know from the previous subsection that lrnin^oo A logPi(Z„ = 2) = 
logE[e~ x ]. So the previous inequality yields 

limsup A logPi(MiiCA„ = |tfnj + 1|Z„ = 2) < limsup A logE[e" s L 5 »J +4£, L««J]. 

Finally, we prove that limsup,^^ A logE[e~ s L<s™j+ 4L L'5™j] < to conclude. Decomposing the 
expectation with < c < E[X] and using 4L^„j — S\& n \ < yields 

E [ e -*L«»j+«-L«»j] < e -<=L*»j +P(4L L5nJ - S lSnl > -c[6n\) 
< e -cL«»J +P(S [5n] <c[8n\). 

As < c < E[A], by standard results of large deviation theory the probability on the right-hand 
side is exponentially small if E[e~ sX ] = E^ -1- *^] < oo for some s > 0. This yields that 

limsup i logPi (MRCA n = [5n\ + 1|Z„ = 2) < 

n— voo 

and thus 

limsup A log Pi (MRCA n > 5n\Z n = 2) < 0. 

n— >-oo 

□ 

7 Examples with two environments : dependence on the 
initial and final population. 

In this section, we focus on the importance of the initial population. 

Example 1 : the limits of A log Pi (Z n = i) and A log Pi (Z,, = j) may be both finite, 
negative but different. Assume that the environment consists of two states qi and q% such 
that 

r:=P(Qi = gi) = l-P(Qi = ga)>0; = 1; q 2 (0) =p, q 2 (2) = 1 -p, 

withp e (0,1). Then 

ilogPi(Z„ = l) = logr, ilogPi(Z n = 2) >max{logr;log[(l-r)2(l-p)p]}. 
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where the term log r comes from a population which stays equal to 1 in the environment sequence 
(qi, q±, qi, • • • ) and the last term comes from a population which stays equal to 2 in the environment 

2(l-p)p \ 

lim ilogP 1 (Z n = l)< lim A logPi(Z„ = 2). 



sequence (q 2 , qi, q 2 , • • • )• Thus if r is chosen small enough (i.e. r < i+2(i-p)p ) 1 



Example 2 : the limits of A logP fc (Z„ = ft) and A logPi(Z„ = ft) with ft > 1 may be both 
finite, negative but different. Actually, in the case without extinction, we have 

lim i IogPi(Z„ = ft) = ^logP(Z 1 = 1), 

n—>oo 

as soon as Pj(3n : Z n = ft) > and the result is immediate. We give here an example with 
ft = 2 and possible extinction. We first observe that such an example is not possible with one 
environment, i.e. in the Galton- Watson case. Then we introduce a simple example in the random 
environment case and check that it is not in contradiction to Theorem 2.1, before considering the 
asymptotic behavior of the probabilities involved. 

Indeed, for Galton- Watson processes with reproduction law q such that q(0) > and q(2) > 0, the 
fact that g(l) > already ensures that the limits of AlogP 2 (Z„ = 2) and A log Pi (Z n = 2) are 
equal. In the case q(l) = 0, we get that 

Pi(Z n+1 = 2) = P x (Zi > 2,Z n+1 = 2) = ^g(z)P,(Z„ = 2) 

i>2 

whereas killing one of the initial individuals starting from 2 and letting the other survive yields : 

¥ 2 (Z n+1 =2) > 2g(0)^<z(*)P 4 (Z„ = 2). 

i>2 

Thus ¥ 2 {Z n+ i = 2) > 2g(0)Pi(Z„ = 2). A converse inequality is clear, so the limits of A logP 2 (Z„ = 
2) and A logPi(Z„ = 2) have to be equal in the Galton- Watson case with possible extinction. 
Thus, we consider two environments to provide an example that the initial population size is also 
of importance even if extinction is possible. More precisely, let the environment consist of the two 
states q\ and q 2 such that 

r:=P(Qi=g 1 ) = l-P(Q 1 = g a )>0 1 

«i(l)=P , qi(a) = i-p, 

9 2 (0)=p , q 2 {2)=p , q 2 (a) = l-2p, 

with p € (0, i) and a > 2. Note ft = 1 ^ Cl(I), so this example doesn't contradict Theorem 2.1 
where it is assumed that the initial population size is in Cl(I). 
To prove that Pi(Z„ = 2) > P 2 (Z„ = 2), we first observe that 

liminf A logP^Zn = 2) > log(rp), (7.1) 

n— >oo 

which comes from a population staying equal to 1 before the last generation in the environment 
sequence (<?i, qi, . . . , qi, q 2 ). 

Next, let us estimate the extinction probability, given the environment. We first observe that any 
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BPVE whose environments are either qi or q 2 is stochastically larger than the Galton- Watson 
process with reproduction law (and unique environment) q 2 . As a consequence, 

¥ x {Z n = Q\£) < ¥ 1 {Z n = 0|Qi = q 2 ...,Q n = q 2 ) < ¥{Z oc = 0|Q a = q 2 , Q 2 = q 2 ,.. .) =: s e a.s. 

It is well-known that s e is given as the first fix point of the generating function f 2 of q 2 : 

Se =fM = p + P S 2 e + (1 - 2p) S a e . 

Let us now estimate s e . For s = 2p, we have 2p > f 2 (2p) = p + 4p 3 + (1 — 2p)2 a p a if a is large 
enough since p < 1/2. Thus s e < 2p if only a is large enough. 
We get then for alH > 1, k < n, 

¥(Z n = 0\£, Z k = i) < 4 < {2 P y a.s. 

Using this estimate and the explicit law of ¥{Z k+ i = • \ Z k — 2, Q k = q%), we obtain a.s. 

¥ 2 {Z n = 2\£,Q k = qi ,Z k = 2) 

= p 2 ¥(Z n = 2\£, Z k+1 =2) + 2(1 - p)p¥(Z n = 2\£, Z k+1 = 1 + a) 
+ (l-p) 2 ¥(Z n = 2\£,Z k+1 =2a) 

= ¥(Z n = 2\£, Z k+1 = 2) (p 2 + 2(1 - p)p{^ + ^P(Z n = 0|£, Z k+1 = 1 + a - 2) 

+ (1 - pf f^V^n = 0|f , ^fe+i = 2a - 2; 
< F(Z n = 2\£, Z k+1 = 2) (p 2 + 2(1 - p)p{^ + ^ + (1 - p) 2 {2p\ 

If p is small enough (depending on a > 2), we get that a.s. 

P 2 (Z„ = 2\£, Q k = qi ,Z k = 2) < ¥{Z n = 2\£, Z k+1 = 2)3p 2 . 
Analogously, if the environment q 2 occurs in generation fc, we get that a.s. 

¥ 2 (Z n = 2\£,Q k = q 2 ,Z k = 2) 

= 2p 2 ¥(Z n = 2\£, Z k+1 - 2) + p 2 ¥(Z n = 2\£, Z k+1 = 4) + 2p{l - 2p)¥{Z n = 2\£, Z k+1 = a) 

+ 2p(l - p)¥{Z n = 2\£, Z k+1 = a + 2) + (1 - 2p) 2 ¥(Z n = 2\£, Z k+1 = 2a) 
<¥(Z n = 2\£,Z k+1 = 2)3p 2 . 

Next, note that the population starting from Zq = 2 is either always > 2 or extinct. Thus in 
each generation, there are at least two individuals and we may apply the estimates above for the 
subtrees emerging in generation k. Finally we get that 

limsup ± logP 2 (Z„ = 2) = limsup ± logE[P 2 (Z fl - 2\£)] < log(3p 2 ). 

n— >-oo n— j-cjo 

We now choose p small enough such that 3p 2 < rp and recall (7.1 ) to get 

limsup i logP 2 (Z„ = 2) < liminf i logPi(Z„ = 2). 

11 — ^00 n 

Finally, we note that this example shows that, as in the the case without extinction in [S], the 
initial population may be of importance for the asymptotic of the probability of staying small, but 



35 



alive. 
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